Dynamics of Fluctuations of Localized Waves 
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We follow the temporal evolution of mesoscopic intensity fluctuations and correlation in strongly 
localized samples. We find an initial burst in relative transmission fluctuations in random one- 
dimensional (ID) samples due to fluctuations in the arrival time of ballistic transmission. Relative 
fluctuations subsequently rise, then drop to a minimum at a time tm, after which they increase 
rapidly in ID simulations and quasi-lD (QID) measurements. For t > 3tm, results in ID and QID 
samples converge towards predictions of a dynamic single parameter scaling model. These results 
reflect the changing number of modes participating appreciably in transmission as the impact of 
longer lived modes grows with time delay. 
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The scaling of the electronic conductance or of the 
transmission of classical waves characterizes the nature 
of transport^'^ and provides a window on the underly- 
ing quasimodes of the random sample'^''*. Such quasi- 
modes, which we will refer to as "modes" , correspond to 
resonances of an open system. An inverse variation of 
transmission with sample thickness indicates that trans- 
port is Ohmic or diffusive. In such samples, modes ex- 
tend throughout the sample and overlap spectrally. In 
contrast, when transmission scales exponentially, modes 
are exponentially localized within the sample and iso- 
lated spectrally. The description of the transport of elec- 
tronic and classical waves are strikingly similar. This 
can be seen in the equivalence expressed in the Lan- 
dauer relation- between the dimensionless conductance 
g= G/{e'^/h), where G is the conductance, e is the elec- 
tron charge, h is Planck's constant, and the optical trans- 
mittance, T, g = T = J2ab -^"h- Here, Tab is the intensity 
or transmission coefficient of the input and output chan- 
nels a and 6, respectively. These channels can be either 
transverse propagation modes external to the sample or 
positions within distinct coherence areas on the sample 
surface with specified polarization. The scaling of g de- 
pends only upon the value of g itself and the dimension- 
ality of the sample. Anderson localization is achieved 
when g < 1—. Localization of classical waves is of special 
interest because waves can be localized purely by wave in- 
terference without the complicating role of the Coulomb 
interaction. In addition, measurements of fiuctuations 
of relative transmission, which are greatly enhanced in 
the localization transition, can be made for single chan- 
nels or for sums over channels in a random ensemble of 
statistically equivalent samples. Indeed, g and the vari- 
ance of fiuctuations of total transmission normalized by 
its ensemble average are inversely related^"—. Thus the 
statistics of steady state transmission, which can be mea- 
sured for classical waves, characterize propagation and 
localization in random media. 

Recently there has been great interest in the dynamics 
in the localization transition. In contrast to measure- 



ments of the scaling of transmission which can track the 
changing impact of weak localization on samples of differ- 
ent sizes, pulsed measurements may provide the changing 
contributions of underlying electromagnetic modes with 
different decay rates in samples of a particular scale. A 
slowdown of the rate of decay of transmission of classical 
waves has been observed in microwavelSil'', optical^® and 
ultrasoundi^ experiments near and beyond the localiza- 
tion threshold. A related reduction has been observed in 
the spreading of matter waves due to Anderson localiza- 
tion of Bose-Einstein condensates in a random potential 
created by ID optical speckle patterns^S. The ensemble 
average of pulsed transmission was measured for local- 
ized ultrasound just beyond the mobility edge in a slab 
of sintered aluminum beads?^ and for strongly localized 
microwave radiation in random dielectric QID samplesii. 
The ultrasound measurements of ensemble averaged in- 
tensity, {I{t)), were well fit by the self-consistent theory 
of localization (SCLT)iSi2ii^ with a renormalized dif- 
fusion coefficient in space and frequency. In strongly 
localized QID samples, however, the decay rate of mi- 
crowave transmission fell below predictions for SCLT— . 
The slowing down of transport at long times was ex- 
plained in terms of a ID dynamic single parameter scal- 
ing (DSPS) model which neglects mode overlap and av- 
erages over the distribution of decay rates and associ- 
ated transmission strengths of localized modes^iii. The 
relative contributions of long-lived localized modes and 
short-lived "necklace states" is central to understanding 
dynamics^"—. Necklace states are formed by the hy- 
bridization of spectrally overlapping localized states to 
form multiply-peaked modes in space which decay rel- 
atively rapidly through the sample boundary. In ad- 
dition to a decay of average transmission, observations 
of a corresponding growth in correlation with time have 
been made in microwave measurements of in diffusive 
samplesSZi but dynamic measurements of mesoscopic phe- 
nomena have not as yet been carried out in localized sam- 
ples. 

Here we report the first studies of the dynamics 
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of mesoscopic fluctuations and correlation for localized 
waves. ID simulations and microwave measurements 
in QID samples reveal a complex temporal variation of 
transmission statistics following an excitation pulse. A 
short jump is observed in relative fluctuations in ID due 
to the variations in the speed of the ballistic wave in 
different configurations. Subsequently, relative fluctua- 
tions rise and then drop to a minimum at a time tm and 
then rise again. The variation of fluctuations is explained 
in terms of the changing effective number of modes that 
contribute to transmission. Fluctuations are enhanced at 
t < tm and t > t„i by the selective contribution to trans- 
mission of short- and long-lived modes, respectively. For 
t > 3tm, results in ID and QID samples converge to- 
wards predictions of the DSPS model for isolated local- 
ized modes. 

Measurements are carried out on QfD samples con- 
tained in 7.3-cm-diameter copper tubes of various 
lengths. The tubes are filled with alumina spheres of 
diameter 0.95 cm and refractive index 3.14 embedded in 
Styrofoam spheres of diameter 1.9 cm at an alumina vol- 
ume filling fraction of 0.068. Measurements are carried 
out in a narrow frequency window from 10-10.24 GHz 
just above the first Mie resonance of the alumina spheres 
in which localization is fostered by near resonant scatter- 
ing and a dip in the density of states^. Spectra of the 
transmitted field are obtained using a vector network an- 
alyzer. Ensembles of sample realizations are created by 
momentarily rotating the sample tube between measure- 
ments. The average intensity localization length in the 
sample is f = 30 cm^''. 

The response to a Gaussian incident field pulse, 
i?o(i) ^ exp{—t^ /2(Tf) exp{i2TTi/ot) is obtained by tak- 
ing the Fourier transform of the product of the spectra 
of the transmitted field and the incident pulse, £{1/) ~ 
exp(— — i/o)^/2cr^), where cr = {2TTat)~^. This gives the 
time-dependent field E{t), and intensity \E{t)\ . Mea- 
surements are made in a number of experimental config- 
urations. In one, transmission of a plane wave produced 
by a horn antenna is detected by a 4-mm wire antenna 
translated on a 2-mm grid over the output surface of 
61-cm-long samples in 200 sample configurations. In an- 
other experimental arrangement, waves are launched and 
detected using conical horns at a distance from the sam- 
ple. 

Cumulant correlation functions of relative intensity, 
Ia{r,t) — Ia{r,t) / {Ia{r,t)) , with displacement on the 
output surface, Ar, C(Ar, t) 5ia{r,t)dia{r+Ar,t) >, 
for a single incident transverse mode, a, are shown in 
Fig. [T]for different time delays, t, in QID samples with 
L ^ 2^. Here 6Ia{r) — Ia{r) — 1 is the deviation of 
relative intensity from its ensemble average of 1 and 
^ = 30 cm^^ is the intensity localization length. The 
correlation function has the same form for localized and 
diffusive waves, C(Ar) = F{Ar)+K{l+F{Ar))^^. Here 
F{Ar, t) = is the square of the normalized field cor- 

relation function and k is the degree of relative intensity 
correlation at displacements for which F = 0, so that. 
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FIG. 1: (Color online) Spatial correlation, C(Ar) for different 
delay times t from the peak of an incident Gaussian pulse in 
a sample of L = 61 cm ~ 2^. C(Ar) falls at early times and 
then rises, while the square of corresponding field correlation 
functions, F{A{r), shown in the inset, are independent of 
time. The dashed vertical lines at Aro = 12 indicate the first 
zero of F at which point, C = k. C(0) = var(/a(r)) and 
C(Aro) = K are shown on the vertical lines on the right and 
left sides of the figure, respectively. 



for example, k — C(Aro) at the first zero of C. In the 
diffusive limit, k — > and C = is correlated over the 
short range of a speckle spot. Fe is the Fourier trans- 
form of the normalized specific intensity, which is the 
angular distribution of transmitted intensity, which does 
not change in time. As a result, F is independent of i, 
as seen in the inset in Fig. [1] k is seen to drop at early 
times before increasing. 

It was not possible to study early time dynamics of cor- 
relation in QID samples with short pulses because of the 
narrow frequency band over which the statistics of prop- 
agation are uniform. However, the response of narrow 
pulses could be calculated in ID simulations. Configura- 
tions of random samples are constructed using the model 
described in Ref. 16. Samples of = L/a layers of 
equal thickness a, are embedded in vacuum with e = 1 
and wave speed c, the speed of light. The dielectric con- 
stant in each layer, Ci, is a uniformly distributed random 
number between 0.3 and 1.7. The intensity localization 
length is ^ = — L/(lnr) = 22a and the central frequency 
is = 0.26c/ail. 

The temporal variation of relative fiuctuations in 
transmission in ID following a Gaussian incident pulse, 
var(/(t)), calculated using spectra of the transmitted 
field just beyond the output surface for an ensemble of 
50,000 configurations is shown in Fig. [H^a). The sharp 
spike in var(/(t)) (Fig. Ufa)) is due to differences in the 
transit time of the leading edge of the ballistic pulse in 
different random configurations (Fig. lUb)), which leads 
to large fluctuations in lit). The time at which the lead- 
ing edge of the pulse arrives at the sample output, tarr, 
determined by the time at which I(t) reaches the value 
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FIG. 2: (Color online) (a) Variance of relative intensity vs. t 
in ID. (b) Initial jump in var(/(t)) is associated with the dif- 
ferent times of arrival tarr of the rising edge oi I(t) in different 
random configurations. The dashed curve is the configuration 
average of I it). The inset shows the arrival time tarr in dif- 
ferent configurations closely matches the calculated ballistic 
time through the sample with elements of thickness a and 
average phase velocity v. 
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FIG. 3: (Color online) Dynamic variance of normalized trans- 
mission for three sample lengths and a fixed pulse width 
obtained from (a) measurements in QID for three sample 
lengths, L = 40 cm ~ 1.33^, L = 50 cm ~ 1.67^ and 
L = Ql cm ~ 2^, and (b) simulations in ID (curves) and 
the results of the DSPS model for t > 2.5tm (open circles). 



10 ^, is seen to closely track the calculated ballistic ar- 
rival time, o,y/ei/c = {N/v){a/c). Once the leading 
edge of the pulse is transmitted, /(i), remains high for a 
time equal to the inverse bandwidth of the exciting pulse 
and relative fluctuations drop to a minimum. After this 
time, the fluctuations from one configuration to another 
have relative maxima at different times as a result of the 
random phasing of contribution to transmission of modes 
of electromagnetic radiation at different frequencies in 
different sample configurations. The variance of fiuctua- 
tions is then inversely proportional to the effective num- 
ber of modes contributing significantly to transmission. 
At first, only the shortest lived modes contribute appre- 
ciably to transmission since longer-lived modes surrender 
their energy slowly. This leads to a peak in var(/(t)) at 
^ 160a/c. After this time, the contribution of longer 
lived modes begins to be felt since energy in the short- 
est lived modes has already leaked from the sample. As 
a result, the effective number of modes contributing to 
transmission increases and var(/(i)) falls. 

The dynamics of fluctuations over a broader time scale 
for measurements in QID and simulations in ID are 
shown in Figs. Measurement in quasi- ID and sim- 
ulation in ID show minima in var(ja(t)) (Fig. EJa)) and 
var(/(t)) (Fig. |31[b)) at a time trm which we find is inde- 
pendent of pulse bandwidth. This indicates that tm re- 



flects a property of modes of the medium. The minimum 
is associated with a transition from a condition in which 
short-lived modes dominate transmission to one in which 
long-lived modes predominate. Presumably, t„i corre- 
sponds to the time within the intermediate time range of 
Fig. |3]at which the largest number of modes participate 
appreciably to transmission. 

Since the rapid increase of var(/a(r, t))) after i„ 
arises from the dominance of long-lived modes, its be- 
havior might be modeled by the DSPS model pro- 
posed in Ref. 16. According to this model, local- 
ized modes peaked in space at depth z into the sam- 
ple have an amplitude A{'^,z) = exp (— 7 |L — 2z|) at 
the output surface, where the Lyapunov exponent, 7 = 
1/2^, is drawn from a Gaussian distribution, P(7) = 
-\/Z727r^exp [— (7 — 7)^/(27/L)]- with a lower cutoff 
limited by the sample length, i.e., 7 > B/L, where B 
is the only adjustable parameter in our model. The time 
response to a Gaussian incident pulse at long times can 
be written as, 



-^DSPs(0 — 
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FIG. 4: (Color online) Comparison of measurements of the 
dynamics of infinite-range correlation in a Q ID sample of L = 
50 cm 1.67^ and calculations of the variance of intensity in 
the DSPS model for a corresponding ID sample. 

where r(7i,Zi) is the decay rate of a localized state lo- 
cated at Zi with a localization length = and Ui is 
the frequency of the localized state. The explicit expres- 
sion for r(7i,Zi) is given in Eq. (3) of Ref. 16. In our 
Monte Carlo simulations of the DSPS model for /DSPs(i) 
, we average over a window of size Aiy — 0.207c/a. We 
assume that the number of states M excited inside the 
window follows the Poisson distribution with a mean 
M equals to AvpL , where the density of states per 
unit length at i^q is p = 1.32c~^, and the frequencies 
i>i (i = 1, 2, • • • , M) are chosen randomly inside the win- 
dow. The results of the DSPS model for t > 2.5tm are 
shown as circles in Fig. EJb). Excellent agreement be- 
tween the DSPS model and ID simulations is found for 
t > 3t„i. Variances of intensity are seen in Fig. [3] to 
be larger for longer samples at early times as is found in 
steady state. But at later times, variances are larger in 
shorter samples. In this case, the density of trajectories 
at a given time is larger than in longer samples and the 
probability of trajectories crossing is greater. 

Though (/) may be compared to either (/^(r)), {la = 
J2r -^".{r)) or (T = J2a'^a) since these are the same 
in ID, the strength of second order statistics in QID 
depends upon the extent of spatial averaging over the 
speckle pattern. The most apt comparison of var(/) 
to second order transmission statistics in QID is to 



var(s = T/(T)), which is equal to the infinite-range cor- 
relator, Koo = {Sia{r)Sia'{r'))^^L^2J^i^i^, where a ^ a', 
and r and r' are two position on the output surface at 
which field correlation vanishes. Short-range correlation 
of the speckle pattern does not contribute to either var(/) 
or Koo- On the one hand, Kqo is independent of the choice 
of input and output transverse mode or position in QID, 
while on the other, there is no transverse intensity vari- 
ation in ID. A DSPS calculation of var(/) and measure- 
ments of Koo in QID for a sample with i/^ ^ 1.67 excited 
by a pulse of width cr = 5 MHz are compared in Fig. 2] 
Good agreement is obtained for t > 700 ns ^ 3t,„. We 
use the same values for the parameters which appear in 
Eq. (3) of Ref. 16, which were obtained there by fitting 
the decay rate at long times. The density of states per 
unit length in this system is p = 8.67 ns/cm. These re- 
sults demonstrate that for long times, measurements in 
QID approach Monte Carlo simulations of the ID DSPS 
theory in corresponding samples. This convergence even 
in samples for which L /£, is not much larger than unity re- 
flects the similar probability distributions for / and s for 
long times. This is in contrast to the log-normal distribu- 
tion of s in steady state predicted only for L/^^ l^^i^° , 
and reflects the dominance of localized modes at long 
times. 

In conclusion, the dynamics of mesoscopic fluctuations 
of localized waves provides a window on the evolving con- 
tributions of short- and long-lived electromagnetic modes 
of the random medium. For t < t„i but somewhat greater 
than iarr, the transmitted energy is due to modes which 
release their energy quickly, while for t > tm & decreas- 
ing subset of long-lived modes contribute substantially to 
transmission leading to increasingly enhanced mesoscopic 
fluctuations. At trm the contributions to transmission of 
necklace states and long-lived localized modes are most 
democratically represented and the variance of fluctua- 
tions is at a minimum. These results show that complex 
mesoscopic transport phenomena for localized waves can 
be clarified by applying a modal analysis. 
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